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We present here the first experimental study of Wigner islands formed by electrons fioating over 
helium. Electrons are trapped electrostatically in a mesoscopic structure covered with a helium film, 
behaving as a quantum dot. By removing electrons one by one, we are able to find the addition 
spectrum, i.e. the energy required to add (or extract) one electron from the trap with occupation 
number A'^. Experimental addition spectra are compared with Monte Carlo simulations for the 
actual trap geometry, confirming the ordered state of electrons over helium in the island. 



I. INTRODUCTION 

Electrons collect on a two-dimensional sheet when they 
are spread over the surface of a liquid helium film. This 
is due to the attractive charge which appears by polari- 
sation of the helium. Due to the weakness of the image 
charge, electrons sit in the vacuum far away from the sur- 
face and move freely at a fixed distance over the film. As 
such, they constitute a very clean and predictable system. 

When temperature is reduced, kinetic energy decreases 
relative to Coulomb energy; correlations begin to domi- 
nate the electronic structure. Wigner predicted in 1934 
that a phase transition would take place in the infinitely 
extended system, leading to the formation of a two- 
dimensional (2D) electron lattice.'^'Wigner crystallisation 
into a triangular lattice of electrons over helium has been 
observed first by Grimes and Adams. l^EI This phase tran- 
sition takes place below a transition temperature that is 
a function of electronic density. For the infinite system in 
the classical limit,Elin which particles can be treated in- 
dividually and which is mostly appropriate for electrons 
over helium, the transition occurs when the ratio F of the 
average Coulomb interaction energy Eq to the thermal 
energy becomes greater than 137. In the quantum 
case, eventually reached at low enough temperatures, the 
Wigner solid is predicted to form when the Brueckner 
parameter Vg, the interelectron distance normalised by 
the Bohr radius, ob, becomes smaller than 37. Recently, 
somewhat more complex phases and ordering transitions 
have been predicted around that rg threshold .EEI 

An assembly of rectilinear vortices provides another 
example of Wigner lattice formation.!^ Abrikosov's tri- 
angular lattice for vortices induced by magnetic fields 
in type II superconductor^ was first observed indirectly 
by neutron diffractioiJ^^ and later visualised directly by 
electron microscopy imaging by Bitter decoration of the 
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trapped flux lines with ferromagnetic microparticles.'iil 

Also, electrons in semiconductor heterostructures lo- 
calised by an applied magnetic field undergo a magnet- 
ically induced Wigner transition.G^ They form a well- 
studied and well-understood system, the quantum dot, 
which has been the object of extensive research for the 
past two decades .1^2] 



A. Confined electron structures 

When confining 2D-electrons to a restricted planar 
area, the breaking of translation invariance brings in 
important changes with respect to the infinite geom- 
etry case, in particular because thermodynamic phase 
transitions are suppressed by fluctuations. This prob- 
lem of 2D ordering in restricted geometries has been 
extensively studied theoretically, mostly for parabolic 
traps with Coulomb interaction between the electrons, 
first in the limit where electrons behave as classical 
particles,I^ElIlIIll] next in the com putationally more de- 
manding quantum case.'^^'l^ElIIS] 

Finite-size effects become prevalent for N < 100. Such 
is the case in the work described below. New features 
appear that depend on the competition between the tri- 
angular lattice, which takes over for sufficiently large sys- 
tems, and the shape and strength of the confining poten- 
tial, which tends to suppress it. For hard confining walls 
and a flat trap bottom, it is predicted that the ordering is 
mostly affected close to the boundaries, electrons in the 
interior of the 2D island retaining the triangular lattice 
structure. For traps with a parabolic confining potential, 
the particles arrange themselves in circular shells with 
widths small compared to the radius of the shell. Thus, 
even with very few electrons, ordered structures can ap- 
pear, the Wigner molecules. 

The direct observation of such structures in re- 
stricted geometry has been achieved only for systems 
of macroscopic charged particles,l2II and for v ortices in 
superconductor^^] and superfluid helium. I^SIMI 

For electrons over helium, finite islands of electrons can 
now be realisecP^so that well-controlled experiments can 
be conducted on this very clean system, which we report 
here. Direct (visual) observation of the electronic struc- 
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ture cannot be performed in the present experiments: the 
geometric arrangement of the electrons in the island must 
be deduced from properties such as their escape energy 
from the island. 



B. Wigner molecules 

Contrarily to solid-state quantum dots, for which in- 
teractions between electrons decrease exponentially due 
to the screening from surrounding electrodes, electrons 
over helium are located far from conducting bodies; the 
Coulomb interaction is mostly unscreened and gives rise 
to strong, long-range, interparticle correlations. Also, 
there are no nearby impurity, no effective mass correc- 
tion; image charges are well-defined. These features con- 
cur to make Wigner islands of electrons over helium an 
ideal model for the study of strongly correlated few-body 
fermionic systems. 

The electron states and energies depend markedly on 
the trap geometry and on the details of the interparti- 
cle interaction, which in turn depends on the number of 
electrons in the trap and their arrangement. For traps 
with steep walls and fiat bottoms, the electrons tend to 
form a triangular lattice in the island interior and adjust 
to the walls at the trap periphery with a disorganised 
layer. For circular parabolic traps, electrons order first 
on concentric shells following a Hund-type law.l^ These 
structures are referred to as Wigner molecules. 

When fluctuations, either thermal or quantum, are fur- 
ther reduced by lowering the temperature or by adjusting 
the electron areal density, the angular positions of the 
various shells become locked: orientational order sets in 
amongst the previously radially ordered electrons. Some 
particularly stable configurations are found that corre- 
spond to "magic numbers" of electrons with high ordering 
temperatures while other trap occupation numbers lead 
to less stable structures and lower ordering temperatures. 



C. "Phase" diagram 

The various ordering processes taking place in 2D- 
clusters of electrons - the Wigner islands - and the for- 
mation of radially correlated structures - the Wigner 
molecules - have been the subject of a large number of 
theoretical studies. Our interest here lies in systems with 
a small number of electrons, N < 20, confined in rela- 
tively large traps so that the electronic density is low. 

The behaviour of N electrons freely suspended in vac- 
uum over liquid helium and confined laterally in a circular 
parabolic trap by a potential ^rricUjQr'^, rric being the bare 
electron mass and ujq/2tt the harmonic trap frequency, is 
described by the following Hamiltonian 

2—1 ^ ^ 1<J 




Figure 1: Radial (RM) and orientational (OM) "melting" in 
Wigner islands sketched in the n vs kBT/Ec plane for small 
electron occupation numbers A'^ - (0) 10 - (•) 11 - (□) 12 - 
(A) 19 - (■) 20, following Filinov et al.SS The OM bound- 
aries for A''=10, 11 and 20 are very close to the origin and 
cannot be shown. The horizontal (vertical) dash-dash lines 
indicate the quantum (classical) Wigner crystallisation in the 
infinite system, for r-g = 1/n^ = 37 and Fao ~ Ec/ksT = 137 
respectively. 

The terms in curly brackets describe noninteracting 
electrons in the trap. The characteristic length Iq — 
{h/m.cUJo)^^'^ associated with these two terms can be 
viewed as the spatial extent of the electron motion. The 
Coulomb interaction, the last term in Eq.([T|, eventually 
localises and orders the electrons within the trap to dis- 
tances of the order of ro such that Ec = e^/iireoro — 
^nicUjQrQ. Following Filinov et al]^, we can therefore 
take the quantity n = v^/g/rQ, which represents the 
fraction of the trap area actually occupied by the elec- 
trons, a dimensionless measure of the electronic areal 
density. 

The configuration and energies of Wigner islands with 
N electrons are described by Eq.([l|, the eigenstates of 
which however can be found analytically only for N < 
2^ Various approximatio n schem es have to be used for 
higher occupation number.t232M9] Larger electron islands 
require the recourse to numerical simulations (see, for a 
review, Reimann and ManninerP^ . Interactions with the 
environment - heat bath and measuring equipment - are 
not included in Eq.([T]|. In the simulations of Filinov et 
al.Il^the temperature comes in as an input in the Monte- 
Carlo simulations. 

A generalised "phase" diagram for the various order- 
ing processes - radial and orientational - taking place 
in Wigner islands with few electrons has been con- 
structed numerically by Filinov et al.li^whose findings are 
sketched in Figjl] These authors characterise the onset 
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of order by looking at the correlations between electrons 
both radially and angular-wise. In the ordered state, they 
find that vq, introduced above on dimensional grounds, 
is quite close to the mean interparticle distance defined 
by the first maximum of the pair distribution function. 
When coming from the dense, hot (large n or/and T) re- 
gions of the diagram, where the 2D electron cluster is in 
a disordered (liquid or gas- like) state and moving to more 
dilute, cooler region, boundary lines are found first to the 
radially ordered state (RM), and, moving further in, to 
complete ordering where the orientation of the various 
shells become frozen with respect to one another (OM). 



D. Electronic configurations 

As seen in FigJTJ Wigner ordering takes place at low 
densities n and low reduced temperatures k^T/Ec along 
boundaries that differ markedly for various occupation 
numbers TV. Some clusters are more stable than others. 
Specifically, clusters with N =11 and 20 in Figjl] melt 
radially at higher T, n, but orientationally at much lower 
values of T, n.^^^The orientationally more stable clusters 
are called magic clusters, 7V=10, 12, 16, 19 in Ref.jlS). 
and possess a fully frozen structure that is closest to the 
Wigner triangular lattice so that the shells cannot rotate 
on one another. Those that can, A''=ll, 20, are far less 
stable orientationally, but a little more stable radially. 

Somewhat surprisingly, these more radially stable con- 
figurations turn out to be even more stable than the ho- 
mogeneous Wigner crystal itself, which melts in the clas- 
sical limit for temperatures above that set by Foo — 137. 
For the full quantum case, the critical melting density 
for the homogeneous system is set by the value of the 
Brueckner parameter rg = ro/ae = 37. The Brueckner 
parameter can be extended to finite systems by taking 
the definition of rg in an electron cluster given above, 
which yields rg — 1/n?. All boundary lines in Figjl] lie 
below the melting curve for the quantum Wigner crystal. 

It should be mentioned however that these boundaries 
do not correspond to sharp transitions between electronic 
structures with different types of order. The localisation 
process of the electrons into organised structures takes 
place gradually only, as stressed by Ghosal et al.^in par- 
ticular. 

Also, considerably lower critical values for r^, down to 
4, corresponding to much higher critical densities, have 
been reported for clusters with few electrons.^i^ However, 
these findings seem to be an artifact of the computational 
scheme and are not confirmed in subsequent work.'^'^filSil 
It thus does appear, as shown in Figjl] that the critical 
density beyond which quantum fiuctuations dominate, 
destroying the ordered phase of 2D confined structures 
remains of the same order of magnitude as for the open- 
geometry situation when the number of particles, N, is 
large. 

Detailed studies of the ordered configuration of these 
few-body clusters have been performed, mostly by nu- 



merical simulations, both for systems of vortices in rotat- 
ing helium and superconducting disks,^ -'^and in quan- 
tu m dots by a number of authors, notably Peeters et 
al! i4 | i5 | 33 ] |.jjg classical limit, to whose work we com- 
pare our own simulations in Sec jIV| below. 

Visual observations of the ordering of small charged 
metallic spheres confined in circular trap on a plane have 
been performed by Saint- Jean et al.l^ Somewhat unex- 
pectedly, their observations of the ground state config- 
uration of clusters with N up to 30 led to significant 
differences with the predictions for some of the struc- 
tures simulated by more sophisticated Monte-Carlo (MC) 
and molecular dynamics (MD) techniques.^'* ^'^ They 
however agreed with the structures found in older work 
based on a spatial step-by-step minimisation of the free 
energy;^ For instance, the cluster with = 15 is found 
to order in the ground state with an inner shell of 4 
electron^ against 5 electrons in MC-MD calculations. In 
other words, an outer shell with A^ = 10 is favoured by 
the Monte Carlo approach compared to one with A^ = 11 
in the spatial minimisation procedure. Likewise, for a 
cluster with A^ = 16, the observed structure has two 
shells, 5 electrons on the inner shell and 11 on the outer. 
The MC-MD structure would comprise three shells, one 
electron at the centre, 5 and 10 on the two other shells. 

The MC approach has been reexamined by Kong 
et al.lSSl in order to resolve this discrepancy. These au- 
thors paid careful attention to the metastable states and 
the saddle paths between them. They confirm the find- 
ings of earlier work with the same computational tech- 
nique (see also Kong et aLlSS). Noting that the energy 
difference between the true ground state for A^ = 9 and 
16 - in their calculation - and the closest metastable state 
- observed as ground state by Saint-Jean et al.l^is quite 
small, these authors suggest that either the experimen- 
tal situation differs from the one described by the sim- 
ulation, or the observed configuration gets consistently 
stuck in the metastable state. The precise cause for the 
discrepancy is not really uncovered yet although there 
are hints that electron structures are quite sensitive to 
the exact shape of the potential: a departure from a 
parabolic confinement potential going as y ~ r" with 
n = 2 to one with n = 2.2 is enough to alter the computed 
structures and restore agreement with experiment .^SfilSl] 
We shall thus consider that MC calculations carried out 
with the actual confining potential do lead to results in 
agreement with experiments. 



Figure 2: SEM picture of the complete micro-fabricated device (left) and magnified view of the electron trap (right). 



E. Addition spectra 

Addition spectra are built from the energy required to 
add one extra electron in the trap. This energy equals the 
difference in chemical potentials between clusters with N 
and iV + 1 electrons: 

A'{N) = ^l{N + l)-^liN) 

= E{N + 1)-2E{N) + E{N~1) , ^' 

E{N) being the ground state energy of the cluster with 
TV electrons. Once this last quantity has been determined 
numerically, the values of the experimentally more acces- 
sible quantity A^(7V) are known. A study of the details 
of this spectrum thus gives a map of the A^-electron clus- 
ter ground state energies and provides clues about the 
occurrence of ordering in the Wigner island. 

As it became appreciated that addition spectra provide 
experimentally accessible signatures of the onset for the 
formation of shells in Wigner molecules pSl detailed theo- 
retical predictions for various confinement potentials ap- 
peared in the literature (see the reviews by Kouvenhoven 
et al.E^ or Reimann and Manninen'^Sl and, for more re- 
cent work, notably Ghosal et aL^^and Giiglii et al.Sll). 

In the weakly-interacting case, when the last term of 
the right hand side of Eq.([l]) gives a negligible contribu- 
tion, the quantity A^(7V) remains small as shells simply 
fill in and ^{N + 1) ~ f^iN). When adding one electron 
gets a new shell started, for = 3, 6, 10, 15, . . ., the 
addition energy shows spikes that reveal the shell struc- 
ture of the 2D harmonic confinement potential. 

As the confinement potential decreases with respect to 
the Coulomb interaction contribution, i.e, as the electron 
density n diminishes, the Wigner triangular lattice struc- 
ture eventually forces its imprint on the electron cluster. 
A new addition energy signature appears with peaks at 
TV = 3, 7, 11, 13, . . ..31' The crossover occurs around 

— 1/2 

n — rs ~ 0.22: denser clusters are in a liquid-like 
state, thinner ones gradually develop a crystal-like state. 
As emphasised by Ghosal et al. inithe transition is smooth; 



addition spectra are expected to change gradually from 
shell formations to full spatial ordering. The "phase" di- 
agram in Fig jl] shows different thresholds for the onset of 
order as revealed by the study of pair correlations.DS It 
however exhibits the same qualitative trend in the order- 
ing sequence of the electron islands, which can be studied 
by measuring the addition spectra. 

Our goal in the work reported here is a comparison be- 
tween the addition spectra observed in the trap for elec- 
trons over helium used experimentally and the outcome 
of numerical simulations that we have performed using 
the known geometry of the trap. This trap is described 
in Secjn] the experimental procedure in Sec |III| the ob- 
tained spectra in Sec jIV[ and the numerical simulations 
in SecEl 



II. THE ELECTRON TRAP 

The trap used in the present work to confine electrons 
to a 2D island is simple compared to other traps used for 
charged particles, which are usually closed by radiofre- 
quency fields. Here, the electron motion is restricted 
along a plane, on the one side, by the free surface of 
liquid helium that presents an energy barrier of ~ 1 eV 
to the electron, on the other side, by the image charge in 
the fluid bulk that provides electrostatic attraction. 

At temperatures below 1 K the saturated vapour pres- 
sure of helium is very low: electrons float above the sur- 
face of helium in near absolute vacuum. Unlike solids, 
liquid helium contains no impurities. Its topological de- 
fects - the vortices - and its elementary excitations - the 
phonons, rotons and ripplons - are well known, allowing 
relatively simple calculations of interacti ons betw een the 
object in the trap and the environment liSJisJMI ^ com- 
prehensive review of the physics of electrons on liquid 
helium can be found in Ref.(45J. 

The free electron is attracted toward the helium surface 
by its image charge and is repelled by the 1 eV barrier for 
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entering the liquid bulk; it forms what can be viewed as a 
one-dimensional hydrogen-like atom. The dielectric con- 
stant of helium is low (e — 1.057) and the image charge 
Q = e(e— l)/(e-|- 1), e being the electron charge, is small. 
The characteristic lengthscale analogous to the Bohr ra- 
dius is Og — 4aB(e+ l)/(e— 1) « 76 A, mo being the bare 
electron mass, is much larger than the atomic scale. The 
electron in its ground state is floating at l.Sog « 114 A 
above the helium surface. An electric field can be applied 
externally to press the electron on the surface and tune 
both its height and energy. 

If unconfined laterally, the electron moves freely over 
the surface of bulk helium. Its mobility is the highest of 
all condensed matter systems.!^ The transverse motion 
of the electron can be restricted by a system of electrodes 
located below the helium surface. The micro-fabricated 
device used in the present work is shown in Fig j2] and is 
described in full detail in Ref.|47]. 

The right panel of Figj2]shows the trap where electrons 
form a Wigner island. It is a three dimensional struc- 
ture micro-machined on a silicon wafer consisting of i) a 
system of electrodes designed to hold the electrons in a 
well-defined position over the liquid helium surface, ii) a 
single electron transistor (SET) on a pyramidal island at 
the centre of the trap used to detect the presence of elec- 
tronic charges and, possibly, the excitation state of the 
cluster. The pyramidal shape was made with a 5-angle 
evaporation procedure based on the well-known shadow 
technique. All electrodes are made of niobium except the 
SET which is made of aluminium. More details are given 
in Ref.gH]. 

The trap and the reservoir, shown in the left panel 
of Figj2] are 600 nm lower than the surrounding guard 
electrode. Both are filled with liquid helium. The long 
rectangular region acts as a reservoir for surface electrons 
storage. If all electrons in the trap happen to become 
lost, the trap can be replenished by tapping the reservoir. 
The guard electrode is made out of a thick (~0.25 pm) 
layer of Nb deposited on an insulating layer of compa- 
rable thickness. This electrode is used to support, by 
surface tension, the helium film over the reservoir and 
the trap ring, so that the liquid depth is ~ 0.5 pm. The 
electrode itself and the rest of the sample are covered 
by a thin (~200-400 A) film of heHum, held by Van der 
Waals attraction. The bottom of the reservoir is covered 
by two electrodes, made out of thin niobium. These sep- 
arate electrodes control the potentials of the right and 
left halves of the reservoir. They are used to shuffle elec- 
trons from one side to the other. The resulting change in 
capacitance provides a mean of gauging the total mobile 
charge contained in the reservoir. 

A gorge through the guard electrodes connects the 
reservoir to the trap. Usually, the guard electrode is bi- 
ased to negative potential relative to the SET. A poten- 
tial barrier thus forms at the gorge, separating the trap 
and the reservoir. As the right reservoir electrode pro- 
trudes inside the gorge, the potential on this electrode 
strongly affects the height of this barrier. Typical po- 
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Figure 3: Profile of the electrostatic potential in the plane 
of symmetry of the trap for Vsbt = 0.5, 0.2 and V from 
bottom to top. Distances along the a;-axis are in microns 
from the edge of the reservoir. The top of the barrier is in the 
gorge. The difference between the top of the barrier and the 
minimum define the confining potential S. 



tential profiles formed by these electrodes are shown in 
FigjS] All electrodes are capacitively coupled to the SET 
island and can be used as gates. 

The geometric radius of the trap is 3 pm, but its ef- 
fective radius is significantly smaller and depends on the 
configuration of the electrostatic potentials on the var- 
ious neighbouring electrodes. The actual potential well 
is neither parabolic nor axisymmetric. Both shape and 
depth change with the voltage applied to the electrodes. 
It becomes shallower when the lability point, at which the 
confinement completely disappears, is approached. To 
set numbers, a trap frequency of wq/Ztt of ~ 40 Ghz rep- 
resents a low estimate when the trap holds few electrons 
only. The extension of the ground state wavefunction 
of non-interacting electrons is then Iq — {h/nicUJo)^^^ ~ 
21.5 nm. The mean distance between interacting elec- 
trons ro fixed by the Coulomb energy and the parabolic 
confining energy is 200 nm, the average density parame- 
ter n, 0.016, the Coulomb energy Eq, 83.5 K. At a tem- 
perature of 0.2 K, = k^T/Ec = 0.0024. Referring to 
the "phase diagram" in Figjl] the electron island lies well 
into the radially ordered phase region, and, moreover, in 
the classical part of that region. 

Such electrostatic traps have already been constructed 
and operated by Glasson et al.SS Building on the experi- 
ence gained in their work, we have improved the design 
in two ways: 1) the size is made smaller; the electron is 
better confined; 2) the central electrode protrudes from 
the bottom with a pyramidal shape in order to increase 
the confinement and the coupling with the electrons (in 
Figj2] right panel), thus enhancing the detection sensi- 
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Figure 4: Reduced charge on the SET island in terms of the 
potential of the reservoir electrode for several sweeps of Kes- 
The SET potential is +0.6 V, the guard electrode potential, 
0. The inset shows the final steps of the staircase for low 
occupation numbers down to zero, corrected for baseline drift. 
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Figure 5: Variation of the SET output signal U when the 
reservoir voltage, Vies, is swept. The Coulomb blockade os- 
cillation is interrupted by a change of the charge in the SET 
island at ~ 0.5 volt as marked by the arrow. 



pressure becomes high enough. A high vapour pressure 
is required both to ignite the discharge and to thermaHse 
the electrons so that their energy is lower than the energy 
barrier of 1 eV needed to penetrate into the liquid. The 
discharge is ignited by applying ~-500 volts to a wire 
terminated in the middle of the discharge chamber. The 
typical discharge current is ^--^0.1-0.2 pA. The presence of 
electrons is easily detected by applying J/cxc ~100 pvolt 
of 100 kHz voltage to the right reservoir electrode and 
measuring the voltage induced on the left electrode with a 
lock-in amplifier. When electrons appear on the surface, 
the signal changes by 1-20 10~® volt rms. 

After the electrons have been generated and have scat- 
tered over the cell a negative potential is applied to the 
guard electrode while a positive potential is applied to the 
reservoir (the potential difference is typically between 
and +1 V during the discharge). These confining po- 
tentials localise the electrons primarily over the reservoir 
and the ring shaped trap. 

At low temperature, the electrons over the thin Van 
der Waals film become localised while electrons over the 
reservoir and the trap ring, floating much further away 
from the solid substrate, remain mobile. The thickness 
of the helium film covering the reservoir and the trap 
depends somewhat on the amount of helium in the cell. 
As the level of helium decreases, the film thickness also 
decreases. When the helium level drops to '--^ 5 mm below 
the sample surface, only a Van der Waals film remains at 
the centre of the reservoir. In a cell with a small volume, 
it is rather difficult to precisely meter the liquid level 
as significant amounts of helium can remain trapped in 
the fill line by surface tension and fountain pressure. To 
alleviate this problem we have increased the volume of 
the cell below the sample. We determine the amount 
required to fill the cell up to the chip level by measuring 
the capacitance ( ~ 1 pF) between two pads on the chip 
as a function of the volume of helium condensed into the 
cell. When helium liquid starts covering the chip, the 
capacitance increases. We then empty the cell and refill 
it with a quantity of gas corresponding to that just before 
the capacitance increase. 



B. SET readout 



tivity. 



III. EXPERIMENTAL PROCEDURE 

A. Electron seeding and monitoring. 

We seed the electrons on the helium surface by igniting 
a corona discharge in a small chamber separated from 
the rest of the cell by a metallic grid with mesh of a 
few tens of microns. Before the discharge, the cell is 
heated to a temperature ~1.1 K at which the vapour 



The single-electron transistor (SET) located at the 
bottom of the trap o perate s as a sensitive electrometer 
and quantum amplifiei'^fiEIlto detect the presence of elec- 
trons in the trap and, possibly, their quantum state. This 
SET is current-biased and polarised close to the Joseph- 
son quasiparticle peak. 

We record the Coulomb blockade oscillations obtained 
by sweeping the reservoir potential. The voltage across 
the SET, Vg, is modulated according to the total charge 
on its island, which reads 

gi-^Ct/. + ^qo,, . (3) 

i 3 



7 



The first sum is over all the conductors in the system 
with capacitance to the island Ci and potential Ui. The 
second sum runs on the charges induced on the island by 
free charges in the system, such as charges or dipoles in 
the substrate and other electrons over helium. 

Voltage Vs((7i) across the SET is a periodic function of 
charge qi with period e, the electron charge. However, it 
also depends on the bias current and the electronic con- 
figuration in the cell, and varies from run to run. We de- 
termine the functional dependence of VsCffl) on qi exper- 
imentally by sweeping the potential [4 of one of the elec- 
trodes, usually the one that is swept in subsequent mea- 
surements. We then select a portion of the sweep during 
which the background charge distribution did not change 
so that several cycles of the function V{q\) = V{Uc) can 
be superimposed by transforming Uc modulo P, where 
P is an appropriately chosen period. After the mod- 
ulo transformation the selected piece of data is averaged 
and interpolated by a smoothing spline function spl {Uc). 
The rest of the data is fit piecewise with the function 
A spl {Uc + qP), where the amplitude A and the phase q 
are fitting parameters. The amplitude has to be taken 
as a free fitting parameter because the experimentally 
observed V{Ue) varies somewhat in amplitude with Ue- 

The phase q is the charge, expressed as a fraction of e, 
induced on the SET island by free charges in the system. 
The expected value is easy to find from the reciprocity 
principle: it is equal to that induced at the location of 
an electron when the SET island is biased with unit po- 
tential and all the other conductors in the system are 
grounded. We have calculated this potential numerically, 
using the actual 3D geometry of the trap inferred from 
SEM photographs, to obtain a value for the charge q of 
0.5 e. 

To reduce the noise due to fluctuations of the dc volt- 
age across the SET we apply a low-frequency (80-150 Hz) 
modulation with amplitude ~ 100 pvolts to the guard 
electrode. The amplitude of the detected signal is pro- 
portional to the derivative of the periodic function Vs{qi). 



C. Trapping electrons 

After electrons are seeded, we let the system cool down 
from 1 K and proceed to load the trap with a given num- 
ber of them. Typically, at this stage, the guard elec- 
trode is biased to a negative potential, between —0.1 and 
—0.5 volt, and the SET is biased to a positive poten- 
tial between and 0.5 volt. First, we fill the trap with 
electrons by lowering the voltage applied to the reservoir 
electrodes. The electrons are repelled towards the trap. 
Care should be taken not to lower the voltage too much, 
since this can cause irreversible loss of electrons. The 
values of the electrode potentials are determined by trial 
and error. Usually, at least one of the reservoir electrodes 
stays more positive than the guard, although we found 
that we can keep the electrons even with both electrodes 
more negative than the guard. Often, the charging poten- 
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Figure 6: SET phase curve variation for the values of the 
potential on the SET Vset of 0.3, 0.5, and 0.8 volt from left 
to right. The dots mark the escape the last electron from the 
trap. 
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Figure 7: Potential on the reservoir for which the last elec- 
tron leaves the trap. The intercept with the a;-axis yields the 
effective contact potential between aluminium and niobium, 
found to be 0.206 ± 0.005 V. 

tials have to be lowered in the course of the experiment 
as the electron reservoir gradually gets depleted. 

We have tried loading the trap in a controlled manner, 
i.e. with a given number of electrons, but found that the 
electrons usually enter in large, unpredictable bunches, 
even when the potentials of the reservoir electrodes are 
swept slowly. This is contrary to the results reported by 
Glasson et al^The reason for this difference is unclear. 
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After the trap is charged, we start sweeping the poten- 
tial of the right reservoir electrode Vi-cs- When the right 
reservoir voltage is low, the barrier is high and the ring 
remains full of electrons. Sweeping up this voltage re- 
duces the barrier height. When the barrier becomes low 
enough, electrons start leaving the ring. This escape sud- 
denly changes the charge on the SET island. The SET 
island charge variation manifests itself as a phase jump 
in the Coulomb blockade oscillations (see Figjs] in which 
the arrow points to the phase jump). 

Referring to the traces in Figj4j the SET detects a 
sudden change in the charge at Vros ~0.48 volt. Since 
the SET measures charge modulo e, we cannot determine 
directly the absolute value of the total charge but we 
know that the trap flooded with electrons has started to 
empty. 

In this initial phase, electrons leave the trap in such 
a way that individual escapes cannot be resolved. But, 
as Vres is raised further and the barrier height decreases, 
clear step-like jumps become visible. Each electron leav- 
ing the trap changes the induced charge on the SET is- 
land. The corresponding jump of the SET island charge, 
AQ, is plotted in Figj4] against l^cs- Its amplitude de- 
pends on the number of the electrons left in the trap 
and amounts to ^ 0.4e when few electrons only are left. 
This value compares well with that obtained from finite- 
element calculations of the electrostatic field using the 
known trap geometry, ~ 0.5e: electrons escape from the 
trap one by one. 

Finally, at 0.93 volt the last electron leaves the trap. 
Between 0.82 and 0.93 volt there is only one electron left 
in the trap. The next step to the left corresponds to two 
electrons in the trap, and so on. 

The variation of T^os between the jumps depends on the 
Coulomb repulsion between the electrons. Stairs length 
and height increase when the number of electrons in the 
trap decrease. Length increase means a larger reduction 
of the barrier to extract one electron (electron-electron 
interactions decrease). Height increase means that the 
leaving electron induces a larger charge on the island, 
i.e., that the leaving electron is closer to the centre of the 
island. 

These results are similar to the results obtained at the 
Royal Holloway in London. Es] The main difference lies in 
significantly better defined steps, which are both "higher" 
and "longer" in the present work. This is due to bet- 
ter coupling with the pyramidal SET, compared to flat- 
island SET used in London, and to the smaller size of the 
trap, leading to stronger repulsion between the electrons. 
The positions of the steps are more stable in our experi- 
ments as well, as seen in the inset of Figj4] and amenable 
to precise quantitative analysis. 



IV. EXPERIMENTAL ADDITION SPECTRA 

In the following, we shift our attention from SET-phase 
variation signals to addition spectra, which are inherently 
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Figure 8: Addition spectrum for Vset=0.5 V (Vset=0.706 
V in the simulation). Black crosses are simulation results. 
Error bars represent the uncertainty on the fit parameters. 
Empty and filled squares, and circles are experimental data 
for different runs and fall nearly on top of one another. 
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Figure 9: Addition spectrum for Vset=0.3 (experimentally) 
(VsET ~ 0.506 V in the simulations). Crosses are simula- 
tion results. Error bars represent the uncertainty on the fit 
parameters. Squares are experimental results 



more reproducible for the following reason. For a given 
set of external parameter values, such as the potentials 
applied to different electrodes, the potential for which 
one electron leaves the trap can be shifted due to mod- 
ifications in the distribution of trapped charges in the 
substrate. This spurious effect is more likely to occur 
shortly after a new corona discharge, which requires a 
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Figure 10: Energy change with number of electrons in a 
parabolic trap: (A), our results based on a discrete parabolic 
profile, plain line, Kong et al.'s resultJ^ (•, vertical scale on 
the right), relative difference between the two calculations. 



"high" temperature (1-1.2 K) and involves a high voltage. 
As addition spectra are the difference between Vrcs{N—l) 
and V,,s{N), namely Afi{N) = V,,s{N)-V,,,{N -1), po- 
tential drifts are removed. 

We note that T^os is modified by the contact poten- 
tial between niobium and aluminium in the cell. This 
potential is determined in the experimental setup as fol- 
lows. We read in Figj6] the value of the right reservoir 
potential for which the last electron leaves the trap for 
different potentials applied to the SET (i.e. different sizes 
of the trap). These values are then plotted against Vset 
in FigjTj The experimental points fall on a straight line 
that cuts the x-axis for V^cs-R = -0.206± 0.005 V. This 
value represents, to a weak correction due to the image 
charge of the remaining electron, the contact potential 
between niobium and aluminium. It must be taken into 
account to obtain the true value of the potential that acts 
on the electrons. 

We now turn to the staircase-like addition spectra 
themselves. We plot the stair length in terms of the num- 
ber of electrons, i.e., the right reservoir potential varia- 
tion necessary to extract one electron. This gives the 
energy profile in terms of electron number. Experimen- 
tally observed addition spectra are given in Figs, [s] and 
|9]for two trap sizes corresponding to Vset=0.5 and 0.3 
volt respectively. 



V. MONTE-CARLO SIMULATIONS 

In order to interpret the details of these experimental 
results, we have carried out Monte Carlo simulations of 
the addition spectra that correspond to the precise shape 
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Figure 11: Computed addition spectra in this work and that 
of Bedanov and Teeters'^, confirmed by Kong et al.'^,(A), 
at the bottom. The relative difference in % between these 
calculations are plotted at the top, (•, vertical scale on the 
right). 



of the confining potential well in the experimental trap. 

These MC simulations are based on the procedure de- 
scribed by Bedanov and Peeters.^Sl In the low density 
limit, which is the regime attained here, the electrons oc- 
cupy a small fraction only of the total area of the trap. 
They are (mostly) distinguishable and can be treated as 
classical particles. This simplification is also made by 
Peeters et al. 

Our trap is not parabolic but its profile can be deter- 
mined by finite elements calculation based on the known 
geometry (see Figjs]). Indoing so, we have attempted to 
reproduce as accurately as possible the pyramidal shape 
of the SET island. We used in the MC simulations val- 
ues of the confining potential profile computed for a finite 
number of nodes whereas Peeters et al. use an analytic 
form. To find the electrostatic energy of an electron at co- 
ordinate Ti in the trap, we fit the potential profile around 
Ti by a 2D parabola between the three closest nodes of 
the calculated profile. 

In order to check our computational procedure, we 
have reproduced Bedanov and Peeters' results for a 
parabolic trap. We discretise the parabolic potential with 
the same number of nodes as with the actual potential for 
our cell. We then perform the same MC simulations to 
find the energy and the configuration of a known number 
of electrons in the trap. Our results, shown in FigjTo] 
fall very well in line with the published results, those of 
Ref.|14) and the more recent results of Kong et al.l^The 
small systematic difference of ~ 0.75 % seen in Fig 10 
may be due to the discretisation of the confinement po- 
tential. We believe that it does not affect the global 
results and that our simulation does find the true ground 
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X [^m] 

Figure 12: Wigner molecule in a parabolic trap with 31 elec- 
trons. Lines represent iso-potential curves, the minimum of 
the parabola is centred on (0,0). Stars represent the electrons. 
The ground state configuration is the configuration (5,11,15) 
as in Kong et al^ 

state configuration. 

To make sure that our simulations are not noisier than 
Peeters et al.'s, we compare addition spectra. Once 
again, our results are in very good agreement with those 
in the literature, as shown in FigjTT] 

Finally, we compare electronic configurations. The 
ground state of the configuration with 31 electrons is one 
of the most difficult to find due to the closeness of the 
first metastable state (the difference is only 0.004%). The 
result, shown in FigjT2| is once again in agreement with 
Peeters et al.'s. We thus are quite confident that our MC 
simulations lead to correct configurations and energies. 

We now turn to the case of addition spectra in the 
actual (non-parabolic) trapping potential, for which we 
need the electron energy. This energy has two parts: an 
electrostatic part and the repulsion due to the Coulomb 
interaction. The electrostatic part can be attractive if 
electrons are in the trap or repulsive if the electrons are 
on the other side of the barrier (see Figjs}. The elec- 
tron energy reads, discarding the kinetic energy term in 
Eq.([lJ: 

N 2^1 

EWVU^-E.VM+ji^^j: (4) 

where V{ri) is the electrostatic energy, obtained by a fi- 
nite elements calculation, the i-th electron coordinate, 
and e the electron charge. The second term represents 
the Coulomb energy between electrons. The bulk part 
of the energy comes from the trapping; only around 10 
% of the total energy in Eq.([4} is due to the Coulomb 
interaction. 

An electron leaves the trap when its energy overcomes 
the confining potential 6, which is the energy difference 
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Figure 13: Evolution of the confining potential 5 in terms 
of the potential applied on the reservoir electrode Kcs- Evo- 
lution of the average energy per electron for 1 electron and 
5 electrons in the dot. Crosses are Monte-Carlo simulation 
results. The plain line is a linear extrapolation of the Monte- 
Carlo results as explained in the text. The crossing point be- 
tween 6 and this average energy gives the potential for which 
an electron leaves the trap. 

between the minimum and the energy at the barrier (see 
Figjs}. That is, a given configuration is stable as long as 
the average energy per electron remains lower than the 
confining potential: 

En = ■ (5) 

This assumes that electrons have identical energies. 

We first compute the evolution of the confining poten- 
tial 6 in terms of the potential applied to the reservoir. 
The contact potential must be taken into account: when 
the potentials applied to the SET electrode and to the 
reservoir electrode are Vset and V^os, the potentials used 
in the simulations are Vset + 0.206 V on the SET elec- 
trode and Vi-cs left unchanged on the reservoir electrode. 
For N electrons in the trap, we compute the total en- 
ergy for different values of the reservoir potential. When 
the barrier too low - i.e. when the reservoir potential 
is too close to the threshold when the electron is about 
to escape - the Monte-Carlo simulations fail to find the 
ground state. The starting temperature in the simula- 
tions is too high and allows electrons to escape readily 
over the barrier. Using a lower starting temperature to 
circumvent the problem does not lead to the real ground 
state of the configuration. To resolve this issue we cal- 
culated the energy for values of the reservoir potential 
slightly lower than the escape threshold, for which the 
simulations did find the ground state, and extrapolated 
the results linearly to higher values of this potential. 

Figure[T3]shows the evolution of the confining potential 
and the average energy per electron for N=l, 3, and 5 in 
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terms of the reservoir voltage. An electron leaves the 
trap when these curves intersect. 

Figures |8] and |9] show the comparison between the out- 
come of these calculations and the experimental results 
for respectively Vset=0.5 and Vset=0.3 V. The size and 
the shape of the trap depend on Vset and the results 
are significantly different. Different potential sweeps fol- 
lowing a given corona discharge are shown on the same 
graph in FigjS] the experiment gives quite reproducible 
results. The scatter on the experimental points is less 
than the uncertainty on the simulations. The observed 
addition spectra are quite reproducible as long as the cell 
is kept cold, below 1 K. When it becomes necessary to 
replenish the electron reservoir, the temperature is raised 
and a new corona discharge ignited. Then the distribu- 
tion of stray charges changes and the addition spectra 
fine structure, which is quite sensitive to the potential 
profile, also changes. 

The addition energy is well predicted for > 3. For 
very few electrons {N ^ 3), the Monte-Carlo simulation 
yields a larger addition energy. 

A recurrent feature of the observed spectra is a peak for 
N = 6 followed by a trough for N = 7. Referring to the 
work of Giiglii et alJ^, this indicates that electrons order 
in shells and not on a triangular lattice, which would give 
a peak at A^ = 7. The corresponding trap frequency is 
~ 60 GHz and the electronic density n ~ 0.017. These 
observations correspond to the phase diagram shown in 
Figjl] according to which, at an electron temperature of 
^ 200mK, the cluster is radially oriented (formation of 
shells) but not orientationally oriented. 

As a rule, the experimental addition energies are found 
smaller than the calculated ones, especially for low occu- 
pation number. This observation probably means that 
the Wigner island is in an excited configuration. Its en- 
ergy is higher than that of the ground state; a lesser de- 
crease of the barrier is required to extract one electron. 
On the contrary, higher values than the calculated ones 
remain unexplained (e.g., A^ = 5 in Figjs}. The extra- 
neous source of noise energy that seems to be present in 
the experiment possibly comes from the tight coupling 
with the SET. We have evidence that the SET backac- 
tion affects the temperature of the electrons in the island 
in a way that depends on the SET bias current, either 



heating or cooling with respect to the bath temperature. 
This effect is under study. 

It is known from the moving pictures of vortex lines 
configuration in a superfiuid rotating bucket by Williams 
and PackardlSSlthat the vortices tend to jump around ran- 
domly. In order to damp the vortex motions so that they 
could be photographed, the authors of Ref.p3] added '^He 
impurities to the *He superfiuid to bring in some dissipa- 
tion. In our Wigner islands, no such dampener is intro- 
duced. Due to the extremely high mobility of electrons 
over helium, it is quite likely that the electronic configu- 
ration is also extremely unsteady. 

VI. CONCLUSION 

We have studied the confinement of a small number 
of electrons in a trap over a liquid helium film. Sta- 
ble configurations down to a single electron can be ob- 
tained reproducibly and their energy recorded with an 
SET readout. The experimental addition spectra com- 
pare well with those obtained in Monte Carlo simula- 
tions in a trapping potential directly derived from the 
cell geometry. Charges trapped in the dielectric parts 
of the sample modify only weakly the potential profile. 
The good agreement between actual experiments and MC 
simulations 1) confirms the validity of the model (and 
assumptions) used in the simulation, 2) shows that no 
uncontrolled, or unforeseen, feature plays a significant 
role in the physical system, 3) opens the way, once the 
problem of repeatability of addition spectra after corona 
discharges is solved, to more detailed studies of these 
Wigner islands, and, in particular, of orientational or- 
dering. 
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